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Abstract 

A highly precise model for the motion of a rigid Earth 
is indispensable to reveal the effects of non-rigidity in 
the rotation of the Earth from observations. To meet 
the accuracy goal of modern theories of Earth rota- 
tion of 1 microarcsecond (pas) it is clear, that for 
such a model also relativistic effects have to be taken 
into account. The largest of these effects is the so 
called geodetic precession. 

In this paper we will describe this effect and the stan- 
dard procedure to deal with it in modeling Earth 
rotation up to now. With our relativistic model of 



Earth rotation ( Klioner et al.l . 120011 ) we are able to 



give a consistent post- Newtonian treatment of the ro- 
tational motion of a rigid Earth in the framework of 
General Relativity. Using this model we show that 
the currently applied standard treatment of geodetic 
precession is not correct. The inconsistency of the 
standard treatment leads to errors in all modern the- 
ories of Earth rotation with a magnitude of up to 200 
fj,as for a time span of one century. 

1 Introduction 

Geodetic precession/nutation is the largest relativis- 
tic effect in Earth rotation. This effect has been dis- 
covered already a fe w years a f ter th e formulation of 
General Relativity ( de Sitter . 19161) and very early 
it was recognized to be important for Earth rota- 
tion. It results mainly in a slow rotation of a geo- 
centric locally inertial reference frame with respect 
to remote celestial objects roughly about the eclip- 



tic normal. Due to its relatively large magnitude of 
about 1.9" per century, which is 3 x 10"'* of the gen- 
eral precession, corresponding corrections are used in 
all standard theories of precession and nutation since 
the lAU 1980 theory (Seidel manni . il982i ). 

The standard way to consider geodetic precession 
in Earth rotation theories up to now was the follow- 
ing: firstly, using purely Newtonian equations, one 
computed the orientation of the Earth in a geocen- 
tric, locally inertial reference frame. To obtain the so- 
lution with respect to the kinematically non-rotating 
Geocentric Celestial Reference System (GCRS) the 
corrections for geodetic precession were then simply 
added, as described for example in iBretagnon et al 
(1997. Section 8). These corrections can be cal 
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culated separately, since they are completely inde- 
pendent of the r otatio nal state of the Earth, e. g. 
Brumberg etall (|l99l[ ). 

The purpose of this work is to demonstrate that 
the standard way of applying the geodetic preces- 
sion is not correct. After stating the problem in the 
following Section, we explain shortly our relativistic 
model of Earth rotation used for this study in Sec- 
tion [31 In Section m we describe how the corrections 
for the geodetic precession can be computed, while 
in Section [S] two different, but equivalent and correct 
ways to obtain a GCRS solution are given. In the 
last section of this paper we compare our solution to 
published ones and draw concluding remarks. 

Throughout the paper we will use the following 
conventions: 

- Lower case Latin indices take the values 1,2,3. 

- Repeated indices imply the Einstein's summa- 
tion irrespective of their positions, e. g. x'y' = 



eabc is the fully antisymmetric Levi-Civita sym- 
bol, defined as Cabc = (a — b){b — c)(c — a)/2. 

Vectors are set boldface and italic: X — X^, 

while matrices are set boldface and upright: P = 

pij 

The choice to use index or vector notation for a 
specific formula is done with regard to readabil- 
ity and clarity. 



where c is the speed of light in the vacuum, G the 
gravitational constant, va the BCRS velocity of the 
body A with mass Ma, ve is the velocity of the geo- 
center, tea the vector from body A to the geocenter 
and rEA its Euclidean norm. The angular velocity 
rJcp corresponds to the orthogonal matrix i?'-' so 
that the respective kinematical Euler equations read 



^GP ~ 2 ^"^^ ^ 



"'iT) — R''%T). (3) 



2 The GCRS and geodetic pre- 
cession 

The Geocentric Celestial Reference System is offi- 
cially adopted by the lAU to be used to describe 
physical phenomena in the vicinity of the Earth and, 
in particular, the rotational motion of the Earth. 
The GCRS is connected with the Barycentric Celes- 
tial Reference System (BCRS) by a generalized ver- 
sion of the Lorentz transformation. This transforma- 
tion was chosen in such a way that the GCRS spa- 
tial coordinates Xk are kinematically non-rotating 
with respect to the BCRS coordinates, i. e. no addi- 
tional spatial rotation of the coordinates is involved 
i n the transformat ion from one system to the other 



( Soffel et al.Ll2003D 



Since the origin of the GCRS coincides with the 
geocenter and the Earth is moving in the gravita- 
tional field of the Solar system a local inertial frame 
with spatial coordinates Xd slowly rotates in the 
GCRS: 



X}, = R^^T)Xi 



K- 



(1) 



Here i?*-' is an orthogonal matrix, the time variable 
T is the Geocentric Coordinate Time TCG. Due to 
this rotation the equations of motion in the GCRS 
contain a Coriolis force. The locally inertial analo- 
gon of the GCRS is called dynamically non-rotating. 
This rotation between the kinematically non-rotating 
GCRS and its dynamically non-rotating counterpart 
is called geodetic precession. The angular velocity of 
geodetic precession rJcp is given by 



1 



GMa 



2 



A ^EA 



■^ve- 2va 



X rEA 



(2) 



This equation can be easily verified by direct substi- 
tution of the matrix elements. 



3 Model of Earth rotation 

A complete and profound discussion of our relativistic 
model of Ear th rotation can be found in lKlioner et al.l 
20101 ). For the purposes of this work we ne- 
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gleet all other relativistic effects except for the geode- 
tic precession. In particular, we neglect relativistic 
torques, relativistic time scales, and relativistic scal- 
ing of various parameters. Then, in dynamically non- 
rotating coordinates X^ the equations of rotational 
motion of the Earth can be written as 

-^SD=Ln, (4) 
Sd = Sd{Pd;A,B,C), (5) 

Ld = Ld{Pd', CirmSim] Xad)- (6) 

Here i_D is the torque and So the angular momen- 
tum in the dynamical non-rotating frame and Pd is a 
time-dependent orthogonal matrix transforming the 
coordinates Xo to a terrestrial reference system Y", 
where the gravitational field of the Earth is constant: 



P3\T)X^. 



(7) 



A,B,C are the principle moments of inertia of the 
Earth, Cim,Sim are the coefficients of the gravita- 
tional field of the Earth in Y and xad are the BCRS 
coordinates xa of body A (Sun, Moon, etc.) rotated 
by the geodetic precession: 



(8) 
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The matrix P/^ can be parametri zed by Euler an- 



les t p, and w in the usual way (jBretagnon et al 



[1997.). Thus, these three angles as functions of time 
T represent a solution of Eqs. (H])-®. Note that the 
only difference between a purely Newtonian solution 
of Earth rotation and Eqs. (H])-® is that the torque 
should be computed by using rotated positions xad 
of external bodies and not the normal BCRS posi- 
tions xa- This reflects the fact that the coordinates 
Xd rotate with respect to the BCRS. 

The corresponding equations in the kinematically 
non-rotating GCRS take the form 



GP 



X S 



d 

dT 

SK^SK{PK;A,B,C;nGp), 

Lk — Lk{Pk', Clm, Slm'jXA)- 



(9) 
(10) 
(11) 



The torque Lk is defined by the same functional form 
as L]j, but the coordinates of external bodies are 
taken directly in the BCRS. Eq. ^ contains an ad- 
ditional Coriolis torque proportional to Oqp . Besides 
this, the angular momentum Sk in the GCRS explic- 
itly depends on the geodeti c precession flgp . Fo r de- 
tails of these equations see Klioner et al. ( 20101 ) and 
references therein. 

From Eqs. ([ij and ([7]) it is clear, that the solutions 
of these two sets of equations are related by 



p ab 



cb 



4 Computing the geodetic pre- 
cession 

To determine the effect of geodetic precession one 
has to compute the matrix R. This can be done by 
a numerical integration of Eqs. ©-(O. It should 
be remarked that care has to be taken how to rep- 
resent this matrix properly. To avoid the disconti- 
nuities that can arise when using the common Euler 
angles to describe an arbitrary rotation, we decided 
to use quaternions to represent this matrix. With 
matrix R, solution and Eq. we can calcu- 

late the differences S(p, Sip and 6uj between the Euler 
angles ip, ip and uj corresponding to matrix V k and 
those corresponding to P/j. 



In the process of computing and verifying the re- 
sults we have found and corrected a sign error for 
the correction induced by the geode tic precession 



for angle tp in iBretagnon et al.l ()1998[ ) . Taking this 



sign error into account the differences between the 
analytical solution for dip, Sij} and 5uj deri v ed b 
Brumberg et al. ( 199l[ ) and IBretagnon et al.l ( 199 
and our numerical solution are below 1 /zas. They 
are shown in Fig. [TJ The remaining differences are 
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Figure 1: Differences (in /ias) for the analyt- 
ical solution fo r geod etic precession derived by 



(12) 

IBretagnon et al. 



tion. 



( 1998[ ) and our numerical solu- 
The sign error in the SMART solution 



( Bretagnon et al. . 1998[) for angle is corrected here. 



explained by the limited accuracy of the analytical 
treatment of this effect by the other authors com- 
pared to our numerical result. 



5 Computing the GCRS solu- 
tion 

According to the equations given in Section [3] there 
are two ways to compute the matrix V k correspond- 
ing to the solution of the rotational motion of the 
Earth with respect to the kinematically non-rotating 
GCRS: 
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1. One can numerically integrate Eqs. dl])-® and 
obtain the solution Fjj with respect to dynamically 
non- rotating coordinates X]j. Then one can correct 
for geodetic precession using the matrix R and Eq. 

to rotate the solution into the GCRS. 

2. One can integrate Eqs. ©-([H]) and directly 
obtain Pk- 

Obviously, the initial conditions for Eqs. (jl])-® and 
([5|)- PT|l are again related by ([T^ taken at the initial 
epoch. 

We have implemented both of the above mentioned 
possibilities to compute Pk and verified that the dif- 
ferences in tp, ip and uj computed in the two ways 
represent only numerical noise at the level of 0.001 
/ias and less after 100 years of integration. 

The implementation is done in an efficient way. 
The numerical integration of the matrix R runs for 
example simultaneously with the numerical integra- 
tion for Pd- The relative running times between 
a purely Newtonian integration and both ways de- 
scribed above are given in TablelH Further details on 
our numerical code and its c a pabili ties can be found 
for example in lKlioner et al. ( 2008[ ). 



200 r 



Table 1: Relative CPU times for various integrations. 

Newtonian case 1.00 
Interation of Pd with rotated ephemeris 1.21 
Direct integration of Pk 1-08 

A purely Newtonian model differs from Eqs. 
only by the positions of the solar system bodies 
used to compute the torque on the Earth: the Newto- 
nian model uses the BCRS ephemeris directly, while 
for Eqs. (jl])-® one has to rotate this ephemeris 
according to Eq. ([8]). It is this rotation that has 
never been considered before in any theory of Earth 
rotation, which represents the main source of incon- 
sistency in the standard way of taking the geodetic 
precession into account. The effect of the rotation 
of the ephemeris on the Euler angles (f, ip and uj is 
shown in Fig. [2j One finds that the error due to this 
inconsistency amounts to 200 /J,as after 100 years of 
integration. 

A summary of the interrelations between the cor- 
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Figure 2: Differences (in /zas) for the Euler angles be- 
tween a purely Newtonian solution and the correct so- 
lution in dynamically non-rotating coordinates. The 
latter is obtained by using ephemeris data rotated 
according to Eq. ([8]). 



rect solutions for the Earth rotation in dynamically 
and kinematically non-rotating coordinates as well 
as Newtonian and "kinematically non-rotating" solu- 
tion derived in the standard, inconsistent way is given 
schematically in Fig. [3] 

6 Difference to existing GCRS 
solutions 

The difference between the Euler angles of the 
GCRS solution obtained in this work and the pub- 
lished kinematicahy no n-rotating SMART solution 
( Bretagnon et al. . 19981 ) is given in Fig.|4l Analysing 
the sources of these differences, one can identify three 
components: 

- influence of the rotation of the ephemeris shown 
in Fig. [2] due to the incorrect treatment of the 
geodetic precession; 

- sign error in the correction for geodetic preces- 
sion in ip; 
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Figure 3: Schematic representation of the differences 
in the standard and the correct way to treat geode- 
tic precession. "GP" stands for geodetic preces- 
sion/nutation. Each gray block represents a solution. 
A solid arrow means: add precomputed geodetic pre- 
cession/nutation to a solution to get a new one. A 
dashed arrow means: recompute a solution with in- 
dicated change in the torque model. 



errors of the analytical SMART solution com- 
pared to the more accurate n umerical integra- 
tion, as already discussed in iBretagnon et al 



( 19981) . 
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Figure 4: Differences (in /las) between the published 
kinematically non-rotating SMART solution and the 
correct kinematically non-rotating solution derived in 
this study. 



ft should be remarked that the above-mentioned in- 
consistency is not only restricted to the SMART so- 
lution, which we used in this study for comparison, 
but is also contained in the lAU 2000A Precession- 
Nutation model as described in section 5.5.1 of 



the lERS Conventions (jMcCarthy and Petitl . 120041 ). 
Therefore Fig. |4] allows us to conclude that the ex- 
isting GCRS solutions for rigid Earth rotation are 
wrong by 1000 /J,as in ip, 200 fias in ijj and 100 /zas in 
uj within 100 years from J2000. It can be shown that 
these differences cannot be eliminated by fitting the 
free parameters of our model, namely the moments of 
inertia, the initial Euler angles and their time deriva- 
tives. 

Let us finally note that the effects of non-rigidity 
in the Earth rotation and the inaccuracies of the cor- 
responding models, e. g. for the atmosphere and 
oceans, are significantly larger than the effects dis- 
cussed in this work. Nevertheless to avoid a wrong 
geophysical interpretation of the observed Earth ori- 
entation parameters, the treatment of geodetic pre- 
cession should be done along the lines presented in 
this paper. 

References 

Bretagnon, P., Francou, G., Rocher, P., and Simon, 
J. L. (1998): SMART97: a new solution for the ro- 
tation of the rigid Earth, Astron. Astrophys. 329, 
329-338. 

Bretagnon, P., Rocher, P., and Simon, J. L. (1997): 
Theory of the rotation of the rigid Earth., Astron. 
Astrophys. 319, 305-317. 

Brumberg, V. A., Bretagnon, P., and Francou, G. 
(1991): Analytical algorithms of relativistic re- 
duction of astronomical observations, in Journees 
1991: Systemes de Reference Spatio-temporels, pp 
141-148. 

de Sitter, W. (1916): Einstein's theory of gravitation 
and its astronomical consequences, Mon. Not. R. 
Astron. Soc. 76, 699-728. 

Klioner, S. A., Gerlach, E., and Soffel, M. H. (2010): 

Relativistic aspects of rotational motion of celestial 



5 



bodies in S. A. Klioncr, P. K. Scidclmann, and 
M. H. Soffel (eds.), lAU Symposium 261, pp 112- 
123. 

Klioner, S. A., Soffel, M. H., and Le Poncin-Lafitte, 
C. (2008): Towards the relativistic theory of pre- 
cession and nutation in Journees 2007: Systemes 
de Reference Spatio-temporels, pp 139-142. 

Klioncr, S. A., Soffel, M. H., Xu, C, and Wu, X. 
(2001): Earth's rotation in the framework of gen- 
eral relativity: rigid multipole m,om,ents in The 
Celestial Reference Frame for the Future (Proc. of 
Journees 2007), N. Capitaine (ed.), Paris Obser- 
vatory, Paris, pp 232-238. 

McCarthy, D. D. and Petit, G. (2004): lERS 
Conventions (2003), lERS Technical Note No.32, 
BKG, Frankfurt. 

Seidelmann, P. K. (1982): 1980 lAU Nutation: The 

Final Report of the I A U Working Group on Nuta- 
tion, Celest. Mech. Dyn. Astron. 27, 79-106. 

Soffel, M. H., Klioner, S. A., Petit, G., Wolf, 
P., Kopeikin, S. M., Bretagnon, P., Brumberg, 
V. A., Capitaine, N., Damour, T., Fukushima, 
T., Guinot, B., Huang, T., Lindegren, L., Ma, 
C, Nordtvedt, K., Ries, J. C, Seidelmann, P. K., 
Vokrouhlicky, D., Will, C. M., and Xu, C. (2003): 
The lAU 2000 Resolutions for Astrometry, Celes- 
tial Mechanics, and Metrology in the Relativistic 
Framework: Explanatory Supplement, Astron. J. 
126, 2687-2706. 



6 



